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Swift/XRT monitoring of the candidate Supergiant Fast X-ray 
Transient IGR J16418-4532 
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ABSTRACT 

We report on the Swift monitoring of the candidate supergiant fast X-ray transient (SFXT) 
IGR J 1641 8— 4532, for which both orbital and spin periods are known (~ 3.7 d and ~ 1250 s, 
respectively). Our observations, for a total of ~ 43 ks, span over three orbital periods and 
represent the most intense and complete sampling of the light curve of this source with a 
sensitive X-ray instrument. With this unique set of observations we can address the nature 
of this transient. By applying the clumpy wind model for blue supergiants to the observed 
X-ray light curve, and assuming a circular orbit, the X-ray emission from this source can be 
explained in terms of the accretion from a spherically symmetric clumpy wind, composed of 
clumps with different masses, ranging from ~ 5 x 10 16 g to 10 21 g. Our data suggest, based 
on the X-ray behaviour, that this is an intermediate SFXT. 

Key words: X-rays: binaries - X-rays: individual (IGR J 1 64 18— 4532) 



> . 1 INTRODUCTION 



" 10 erg cm 2 s 1 



and 4.8 x 10~ 10 erg cnT 2 



The X-ray transient IGR J16418— 4532 was discovered by IN- 
TEGRAL during observations of th e black hole X-ray transient 
4U 1630-47 on 2003 February 1-5 dTomsick et ai1l2004|). wi th a 
20-40 keV flux of 3 x 10~ n erg cm" 2 s~ MSguera et alJfeOOrjl) re- 
ported on further INTEGRAL observations, taken on 2004 February 
26, during which fast (~ 1 hr) X-ray outbursts were observed that 
peaked at ~ 80 mCrab (20-30 keV). This behavior led Sgue ra et alj 
d2006l) to propose that IGR J16418— 4532 i s a candidate super giant 
fast X-ray transient (SFXT). SFXTs (e.g. ISguera et alj|2005l) are 
high mass X-ray binaries (HMXBs) associated with OB supergiant 
stars via optical spectroscopy that show an X-ray dynamic range of 
3-5 orders of magnitude. Indeed, they are characterized by bright 
(peak luminosities of 1 36 -10 37 erg s" 1 ), s hort (a few hours, as ob- 
served by J7Vr£Gft4L. ISguera et alj|2005t iNegueruela et alj|2006t) 
X-ray outbursts significantly shorter than those of typical Be/X- 
ray binaries, and a quiescent lum inosity of ~ 10 32 erg s _1 (e.g. 
lint Zandl2005l ; lBozzo etai]|2010h . 

In their syste matic analysis of th e INTEGRAL observations 
from 2003 to 2009, iDucci et al.1 < l2010h detected 23 outbursts, for 
an activity duty cycle of ~ 1 % (one of the highest among the 14 
SFXTs and SFXT candidates they examined) and fluxes ranging 



between 1.3 x 10 
s" 1 (18-100 keV). 

XMM-Newton observations dWalter et al.1 |2006h showed a 
heavily absorbed (7Y H = (1.0±0.1) x 10 23 cm" 2 ) pulsar (P spin = 
1246± 100 s) with an unabsorbed flux of 1.3 x 10~ 10 erg cm -2 s _1 
(2-100 keV) in a fit performed in conjunction with the INTEGRAL 
data. 

IChatv etail 12008) proposed, based on the XMM-Newton po- 
sition, 2MASS J16415078-4532253 as the best NIR counterpart 
candidate; they also fitted the spectral energy distribution (SED), 
obtaining a temperature of 32 800 K for the massive compan- 
ion, consistent with an OB spectral classification. SED fitting by 
iRahoui et al] d2008h yielded A v = 14.5 mag, and a 08.5 spectral 
type; depending on the luminosity class they derive a distance of 
4.9 (V), 8.3 (III), or 13 kpc (I). 

A periodicity of ~ 3.75 d was discovered b y ICorbet et al 
d2006h in the Swift/Burst Alert Telescope (BAT, Barthel mv et al 



120051) (3.753 ± 0.004 d, MID 53360-53630) and RXTEI ASM 
(3.7389 ± 0.0004 d, MID 50091-53810) light curves, which was 
interpreted as the orbital period of a binary system. The hint of a to- 
tal eclipse at phase 0.5 o bserved in the BAT data is consistent with 
a supergiant companion. iLevine et alj|2010i) refined the orbital pe- 
riod value to 3.73886±0.00028 d, by using the RXTEI 'ASM data up 
to MID 55224. The observed values of the orbital and pulse periods 
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Table 1. Summary of the Swift/XKT observations. 



ObsID 



Start time (UT) 



End time (UT) 



Exp. 

(s) 



0003 1929001 b 

00031929002 

000431050()l c 

00031929003 

00031929004 

00031929005 d 

00031929005 

00031929006 d 

00031929006 

00032037001 

00031929007 

00032037002 

00032037003 

00032037004 

00032037005 

00032037006 

00032037007 



2011-02- 
2011-03- 
2011-06- 
2011-07- 
2011-07- 
2011-07- 
2011-07 
2011-07 
2011-07- 
2011-07- 
2011-07- 
2011-07 
2011-07 
2011-07 
2011-07- 
2011-07- 
2011-07- 



18 01:59:43 
04 00:03:19 
27 05:14:54 

13 09:46:55 

14 12:39:39 

15 06:42:44 

15 06:42:50 

16 10:15:02 
16 10:15:05 

18 21:28:28 

19 01:52:41 

20 05:18:42 

21 10:06:18 

22 15:19:43 

23 11:48:19 

29 12:13:30 

30 01:01:30 



2011-02- 
2011-03- 
2011-06- 
2011-07- 
2011-07- 
2011-07- 
2011-07- 
2011-07- 
2011-07- 
2011-07- 
2011-07- 
2011-07- 
2011-07- 
2011-07- 
2011-07- 
2011-07- 
2011-07- 



18 03:43:26 
04 05:01:57 
27 05:23:56 

13 14:43:57 

14 19:31:56 

15 23:00:58 

15 23:01:57 

16 23:05:36 
16 23:06:58 

19 00:51:31 

19 23:19:57 

20 15:22:58 

21 15:27:54 

22 18:44:57 

23 15:23:56 

29 23:34:57 

30 09:27:56 



1951 

2001 

537 

4756 

4832 

273 

1078 

458 

379 

1947 

2673 

3483 

3761 

2812 

4127 

4650 

4287 



0.27 
0.01 
0.80 
0.16 
0.47 
0.72 
0.72 
0.01 
0.01 
0.62 
0.77 
0.01 
0.31 
0.62 
0.85 
0.50 
0.63 



a Mean phase referred to T cpoch = MJD 53560.2 iLevine et alj20ld) . 
b Corresponding to MJD 55610.08. 
c Serendipitous observation. 
d WT mode data. 



indeed place I GR J16418— 4532 in a region of the Corbet diagram 
dCorbel 19861) where wind-accreting HMXBs generally lie. 

In this paper we analyze all the Swift data collected on 
IGR J16418 — 4532. The first set of observations comprises two 2ks 
target of opportunity (ToO) observations performed on 201 1 Febru- 
ary 18 and March 4, following a MAXI alert of a possible detec- 
tion of IGR J 1 64 18 — 45 32 on 2011 February 17 (preliminary re- 
sults were described in Roma no"et alj|201 Tbh . The second set is 
drawn from a monitoring campaign in 2011 July that spans over 
three orbital periods. This Swift monitoring represents the most in- 
tense and complete sampling along the orbital period of the light 
curve of this source with a sensitive X-ray instrument. 



2 OBSERVATIONS AND DATA REDUCTION 

Table Q] reports th e log of the Swift/X-ray Telescope (XRT, 
iBurrows et al1l2005h observations used in this paper. The Swift ob- 
servations of IGR J16418-4532 during 2011 Febru ary and March 
(00031 929001 and 00031929002) were obtained dRomano et all 
2011b) as two ToO observations following a MAXI alert 
(5609524433, later retracted, as due to GX 340+0) of a possible 
detection of IGR J 1 64 1 8 —4532 on 20 1 1 February 17. The data col- 
lected in 201 1 July were obtained as a ToO monitoring program of 
9 scheduled daily observations, each 5ks long, starting on 2011 
July 13. This strategy was devised in order to cover the light curve 
in different orbital phases while maintaining the observing time per 
day reasonably short, thus not hindering gamma-ray burst (GRB) 
observations. However, most observations were cut short by the 
occurrence of GRBs and the discovery of a new soft gamma-ray 
repeater (starting from obsID 00031929005, see Table HJ. Obser- 
vation 00032037001 was performed to determine the best offset 
from the intended target in order to alleviate the problem of single- 
reflection rings from the nearby source GX 340+0 (these rings 



marginally affect the source while in the faint states, so that the 
analysis of the earlier segments was still feasible, with an appro- 
priate choice of extraction regions) A new pointing (hence a new 
obsID) was adopted afterwards. In our work we also include ob- 
servation 00043105001 that serendipitously covers the region of 
IGRJ16418-4532. The 2011 July campaign lasted 18 days di- 
vided in 12 observations for a total on-source exposure of ~ 39 
ks. 

For our binary orbital period search and analysis we retrieved 
the BAT 'orbit-by-orbit' light curves (averaged over Swift's orbital 
period of ~ 90 m, 15-50 keV) covering the data range from 2005 
February 12 to 2011 July 12 (MJD range 53413-55 754) from the 
BAT Transient Monitor iKrimm et afll2006l . [2008T> pageQ- These 
data were further screened to exclude bad quality points (quality 
flag 1 and 2) and referred to the solar system barycentre (SSB) by 
using the task EARTH2SUN. 

The XRT data were processed with standard procedures 
(XRTPIPELINE vO.12.6), filtering and screening criteria by using 
FTOOLS. Both WT and PC events were considered. The selec- 
tion of eve nt grades was 0-2 a nd 0-12, for WT and PC data, re- 
spectively (Burr ows et alj|2005T) . Source events were accumulated 
within a annular/circular region (depending on whether pile-up cor- 
rection was required or not, respectively) with an outer radius of 20 
pixels (1 pixel ~ 2.36"); and background events were accumulated 
from source-free regions far from the contamination of the single- 
reflection rings from GX 340.0+0. For our timing analysis we con- 
verted the event arrival times to the Solar system barycentre with 
the task BARYCORR. Light curves were created for several values 
of signal-to-noise ratio (SNR) and number of counts per bin; all 
were corrected for PSF losses, vignetting and background. For our 
spectral analysis, we extracted events in the same regions as those 
adopted for the light curve creation; ancillary response files were 
generated with XRTMKARF, to account for different extraction re- 
gions, vignetting, and PSF corrections. We used the latest spectral 
redistribution matrices in CALDB (20110 705). 

The UV/Optical Telescope (UVOT, iRoming et ai]|2005h ob- 
served IGR J16418— 4532 simultaneously with the XRT with the 
'Filter of the Day' (FoD), i.e. the filter chosen for all observations 
to be carried out during a specific day in order to minimize the fil- 
ter wheel usage. The exception is observation 00032037007, which 
was requested in the u filter, as opposed to FoD, to attempt to ob- 
tain a detection in that band. The data analysis was performed using 
the uvotimsum and UVOTSOURCE tasks included in the FTOOLS. 
The latter task calculates the magnitude through aperture photom- 
etry within a circular region and applies specific corrections due to 
the detector characteristics. 

All quoted uncertainties are given at 90 % confidence level 
for one interesting parameter unless otherwise stated. The spectral 
indices are parameterized as F u or v~ a , where F u (erg cm - s~ 
Hz - ) is the flux density as a function of frequency v\ we adopt 
r = a + 1 as the photon index, N(E) oc E~ F (ph cm~ 2 s" 1 
keV -1 ). 



3 RESULTS 

3.1 Orbital period from BAT data 

For the orbital modulation, we considered the BAT data in the 
time range MJD 53414-55754. With a standard folding analy- 

1 http://swift.gsfc.nasa.gov/docs/swift/results/transients/ 



IGR J16418-4532 with Swift 3 





BAT 



0.5 



1 

Phase 



1.5 



Figure 1. BAT light curve folded at P ort> 
MJD 53560.20000. 



3.73886 d and T cpoch 



sis of the Solar system barycentred light curves, we measured 
an orbital period of 3.740 ± 0.002 days ( 1-cr error), which i s 
fully consistent with the value obtained by iLevine et alj bold) . 
P orb = 3.73886 ± 0.0002 8 d. Therefore , herea fter we shall adopt 
the more precise measure bv lLevineetailifeOlOl). and we shall refer 
to r c p 0ch = MJD 53560.20000 dCorbet et alj|2006t) . Fig. Q] shows 
the BAT light curve folded over this period. 

The BAT data show an eclipse centred at <j> ~ 0.55. The 
depth is consistent with a total eclipse, and the eclipse FWHM 
duration is 0.17 ± 0.05 of the orbital period. The folded light 
curve is quite flat; the peak-to-trough pulsed fraction, PF pt = 
(T max - i ? m in)/(i ? max + F mln ), where F max and F mi n are the 
observed background-subtracted count rates at the peak and at the 
minimum, is 78 ± 12 % (1-cr error). 



3.2 X-ray position 

We used 20.4 ks of PC mode data and simultaneous Swift/UVOT 
images to obta i n an astrome trically-corrected position (see 
lEvans et al.|[2009l ; lGoad et~aT1l2007l) of: RA(J2000) = 250.46104, 
Dec(J2000) = -45.54092, which is equivalent to: RA(J2000) = 
16 h 41 m 50^65, Dec(J2000) = -45° 32' 27'.' 3, with an uncer- 
tainty of 1'.'9 (90% confidence levelfj This position is 2'.' 5 from 
the source 2MA SS J16415078-453 2253 and 4'.'3 from the XMM- 
Newton position I Walte r etal.l2006l) , and it allows us to confirm the 
2MASS source as the optical counterpart of IGR J16418-4532. 



3.3 UVOT optical counterpart 

At the XRT position of IGR J16418— 4532, no detection was 
achieved in any UV band down to a limit of wl > 21.78, 
m2 > 21.66, w2 > 21.89 mag, respectively. A marginal de- 
tection (2.9cr) was obtained in the u band, at 22.28 ± 0.40 mag. 
T hese magnitudes ar e on the UVOT photometric system described 
in lPoole et al.1 d2008h . and are not corrected for Galactic or intrinsic 
extinction. At this pos ition, the column due to Galactic absorption 



is 1.59 x 10 22 cm"" ( Kal berla et"ai1l2 005). By adopting the con- 
version into Av dPredehl & SchmitJ 19951) , iV H = 1.79 x IQ 21 A v 



2 This new UVOT-enhanced XRT position takes advantage of the n ew cal- 
ibration of the XRT to UVOT detector co-ordinates of 2011 August tevansl 
1201 ll) ; see, also: |http://www.sw"i ft. ac.uk/reprocessed.php 
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Figure 2. Swift/XRT 0.3-10keV light curve of IGR J16418— 4532 dur- 
ing observation 00031929001, background- subtracted, corrected for pile- 
up, PSF losses, and vignetting (SNR = 5). Time is referred to the start of 
the observation. 



cm 2 , this corresponds to (Cardelli et al. 198^) A u 
the UVOT u band. 



14.2 mag in 



3.4 Spin period search 

We looke d for evidence of t he 1246 ± 100 s spin periodicity re- 
ported by |Waltei~ at al. (2006) in the Solar system barycentred XRT 
light curves. We only considered observations where the source 
was significantly detected (thus excluding obsIDs 00032037001 
and 00032037004). 

We searched each single observation longer than one spin pe- 
riod and the total combined dataset within ±3 00 s (corresponding 
to a ±3cr error reported bv lWalter etall l2006) by applying both a 
Fourier transform and a Z 2 analysis (a method that does not re- 
quire the binning of the data: lBuccheri et aljl983l) . No statistically 
significant signal was detected either in the individual observations 
or in the combined dataset. With respect to the latter, the upper 
limit on the pulsed fraction, (defined as the semi-amplitude of the 
sinusoidal modulation divided by the mean count rate), computed 
according to lVaughan et "alll ll994l) . is 42 % at the 99 % confidence 
level and for periods between 946 and 1546 s. 

This upper limit is only marginally comp atible with the (64 ± 
10) % value measured by [Walter et alj d2006t) . We note however 
that they were able to detect pulsations only during a fraction of 
their XMM-Newton data and in particular while the source was not 
flaring, and that our data were collected in a variety of states of 
brightness and across almost all orbital phases. So our upper limit 
and their measure may not be straightforward to compare. 



3.5 Light curves 

Fig.[2]shows the 0.3-10 keV light curve of IGR J16418-4532 dur- 
ing the 2001 February 18 observation, binned in order to achieve 
a SNR of 5. The source reached ~ 5 counts s _1 during the ini- 
tial flare and then decayed through a series of smaller flares. The 
lowest point is at 0.16 counts s - . When binned at a SNR of 3 the 
dynamical range, during this observation, is of order 50. 

Fig.[3]shows the 0.3-10 keV light curve of IGR J16418-4532 
of the whole 2011 July campaign after bin-by-bin background- 
subtraction and corrections for pile-up, PSF losses, and vignetting, 
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Figure 3. SwiftfKKY 0.3-10 keV light curve of IGR J16418— 4532 during the 201 1 July monitoring program, background-subtracted and corrected for pile-up, 
PSF losses, and vignetting (SNR= 3). Downward-pointing arrows are 3cr upper limits. Different colours (see the electronic edition) mark different observations 



(see Table[T}. The top axis reports the phase (with P or b = 3.73886 d and T t 



epoch 



: MJD 53560.20000). 
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Figure 4. Swift/XSJ 0.3-10keV flux light curve of IGR J16418— 4532, folded at P orb = 3.73886d and T cpoch = MJD 53560.20000. The data were 
collected in 2011 February and June (dark grey filled circles), and 2011 July (black filled circles). Downward-pointing arrows are 3<r upper limits. The open 
blue circle is a detection obtained by combining three observations (00032037001, 00032037004, and snapshots 1-3 of 00032037007). Also included are the 
2004 August XMM data (light grey curve). 



binned at SNR = 3. The light curve starts at phase 0.16 (assum- 
ing a period of P orD = 3.73886 d and an initial epoch T epoc h = 
MJD 53560.20000) and monitors the source state through over 
three full periods. Superimposed on the long-term orbital modu- 
lation, which follows that seen in the BAT data, flaring is observed 
on short time scales, as sh own in Fig. [2] This behaviour has been 
observed in m ost SFXTs Jsidoli et al j|2008l : lRomano et alj|2009bl 
l2011al , l20ld) . 

The lowest point (as observed by XRT) was collected dur- 
ing the eclipse, by combining three observations (00032037001, 



00032037004, and snapshot^ 1-3 of 00032037007, see open blue 
point in Fig.[4j at 2.5 x 10 -3 counts s~\ This corresponds to an un- 
absorbed 0.3-10 keV flux of ~ 6.0 x 10~ 13 ergcnT 2 s~\andtoa 
luminosity of 1.2xl0 34 erg s" 1 (assuming the optical counterpart 
distance of 13kpc). The lowest point during the campaign outside 
the eclipse was recorded on MJD 55763.43 at 0.022 counts s _1 , 



3 Part of an observation with a continuous viewing of the target, with typi- 
cal timescales of 10 3 s. 



IGR J16418-4532 with Swift 5 




Table 2. XRT spectroscopy of IGR J16418-4532. 



0.1 1 
Count rate (0.3-10 keV) 




0.5 1 

Phase 

Figure 5. Swift/XKF hardness ratio (red open circles, 30 counts per bin) 
as a function of count rate (top) and phase (bottom). The black points are 
weighted means. 



corresponding to an unabsorbed 0.3-10 keV flux of ~ 4.7 x 10~ 12 
erg cm -2 s" 1 , and to a luminosity of 9.5xl0 34 erg s" 1 . As the 
peak count rate is reached on MJD 55610.08 at > 8 counts s" 1 
the observed dynamical range of this source is at least 370, (1400 
considering the points within the eclipse). 

Fig-0shows the 0.3-10 keV light curve, folded using P or b = 
3.73886 d and T cpoch = MJD 53560.20000. We adopted count rate 
to flux conversion factors derived from spectral fits of each obser- 
vation (Section [3. 6t and subsequent interpolation across observa- 
tions where fitting was not feasible. We also include the ~ 10 ks 
2004 August 19 XMM-Newton/EPlC-pn observation dWalter et all 
2006), binned at 256 s, and converted into flux by adopting a 
single conversion factor obtained from fitting the mean spectrum 
with an absorbed power-law model (r = 1.40 ± 0.05, A*h = 
(10.5 ± 0.5) x 10 22 cm" 2 , 1-cr errors; x 2 =0.99/169 dof). The 
eclipse lasts A(j> ~ 0.2 and is centred at <f> ~ 0.55, as observed in 
the BAT data (Fig.[Q. 

Fig.|5]shows the hardness ratio (HR defined as H/S, where H 
is the count rate in the 4-10 keV band and S is the count rate in the 
0.3-4 keV band), as a function of count rate (top), and orbital phase 
(bottom), where the light curve was binned to have 30 counts per 
bin (red open circles). A fit to a constant value yields a hardness 
ratio value of 1.01 ± 0.03 and xt = 0.85 for 130 degrees of free- 
dom (dof), so that the HR is shown to be roughly constant over two 
orders of magnitude in count rate. Furthermore, the HR were also 



ObsID 



JVh 
(10 22 cm- 



Flux" 



C-stat (dof) c 



00031929001 
00031929002 
00043105001 
00031929003 
00031929004 
00031929007 
00032037002 
00032037003 
00032037005 



6.77^ 



7 +2.32 
-1.78 

-1.24 
+3.47 
-2.54 
1-1.00 
-0.82 

5.16+ 417 



5.28 
5.20 



3.93 
4.59 
5.65 
5.14 



2.81 
1.24 
0.97 
1.23 
1.00 
+2.80 
1.96 
+ 1.72 
1.39 



a+0.61 
-0.54 
1+0.51 
-0.45 
-+0.99 
-0.87 
-+0.27 
-0.24 

0.09+ ' 87 



1.38: 

1.04 4 
1.35j 
1.35 4 



1.04 
1.05 
0.62 
1.01 



0.75 
+0.40 
0.36 
+0.38 
0.34 
+0.62 
0.55 
+0.45 
0.40 



21.3612 

11.58+^ 
3.85+^- 
r +o. 
-0 

;+0. 

-0 
)+0 

)+0 

-0 

+0 

-0 
,+0. 



8. 27Z 
1.05j 
7.62j 
6.32j 
2.51 4 
3.88"* 



14 

44 
3!) 

.30 

16 

■S3 

65 

58 

21 

17 

72 

09 

56 

54 

32 

31 

4ii 

37 



0.532/18 
0.845/17 
70.05(76) 
0.935/58 
169.46(141) 
0.647/28 
0.831/31 
1.114/16 
0.838/27 



a Unabsorbed 2-10keV fluxes (10" 11 erg cm -2 ! 
b 2— lOkeV luminosities in units of 10 35 erg s — 1 , at 13 kpc. 
c Cash statistics (C-stat). 



weighted across several phase intervals (black points); we cannot 
find evidence of variations of the HR as a function of phase. 



3.6 Spectral analysis 

We first extracted PC spectra for each segment in which a detection 
was obtained and a minimum of ~ 100 source counts were avail- 
able. We did not consider the WT data for our spectroscopic study 
as they were statistically poor and were heavily contaminated by 
the background from the single-reflection rings from GX 340.0+0. 
The data were rebinned with a minimum of 20 counts per energy 
bin to allow y 2 fitting, with the exc eption of low-count-statistics 
spectra, in which Cash dCasrJll979t> statistics and spectra binned 
to 1 count per bin were used, instead. The spectra were fit in the 
0.3-10 keV energy range with a single absorbed power law model. 
More complex models were not required by the data. The results 
are reported in Table|2] We note that the spectral parameters do not 
show significant variability within the large uncertainties, despite 
the large observed variations in flux throughout the campaign. 

Observation 00031929001 (shown as the dark grey filled cir- 
cles at phase 0.27 in Fig. [4] and detailed in Fig. |2j, in particular, 
deserves special attention. Although it was performed as a follow- 
up of a false trigger (later attributed to the nearby GX 340+0), it 
nonetheless found the source in a very active state, the brightest of 
the whole campaign, with its average 2-10 keV unabsorbed flux of 
~ 2xl0" 10 erg cm" 2 s" 1 , corresponding to ~ 4xl0 36 ergs -1 at 
13 kpc). The peak flux was ~ 2 x 10" 9 erg cm" 2 s _1 (~4x 10 37 
ergs -1 ). An absorbed power-law model yielded a high absorbing 
column Nh = (6.8t 2 j) x 10 22 cm" 2 (in excess of the Galac- 
tic one, 1.59 x 10 22 cm" 2 ) and a photon index V = I^IqJ (see 
TableHJ. 



4 DISCUSSION 

IGR J1641 8— 4532 is a candid ate SFXT, based on the hard X-ray 
behaviour dSguera et alj|2006l). for whic h we know both the or- 
bital (Porbit ~ 3. 7 d, ICorbet et al.ll200o1) and the spin (~ 1250 s, 
IWalter et al.ll2006l) periods. These valu es place IGR J16418— 4532 
in the region of the Corbet diagram (Corbet! [19861) of the wind- 
accreting supergiant high-mass X-ray binaries. In this Paper we 
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have taken advantage of the short orbital period which makes 
IGR J16418— 4532 a very good source to monitor with the XRT, 
to study the X-ray properties of this source as a function of the 
orbital phase, in particular, to address the issue of its SFXT nature. 

The Swifi/BAT and XRT folded light curves of 
IGR J16418— 4532 (Fig. □ and Fig. EJ show the presence of 
a dip that is consistent with zero intensity. When interpreted as 
an eclipse, it can be used to infer the nature of the primary of 
the binary system. The BAT lightcurve of IGR J16418— 4532 
(excluding the eclipse) does not clearly show a flux modulation 
with the orbital phase, suggesting that the eccentricity is not 
large. This is in agreement with the observation that in many 
X-ray binary systems with nearly circular orbits and short orbital 
periods (^4 days), the eccentric orbit produced by the explosion 
of the supernova is circularized on timescale of 10 B yr (consistent 
with the evolution ti mescale of an O star), due to viscous tidal 
interactions (see e.g. IPress et al.|[l975l) . Thus, for simplicity, we 
assume a circular orbit (P or b = 3.73886 d). Let us adopt the 
typical values of stellar masses and radii dMartins et al. I [20051) 
for a primary of spectral type 08.5 and luminosity classes I, III, 
or V, and for the neutron star a mass of Mns = 1.4 Mq, and a 
radius of i?NS = 10 k m. By combining Keple r's third law with the 
relationship given by Rappaport & Josj dl983l) for the radius of the 
primary R, 



(R/a) 2 = cos 2 i + sin 2 i sin 2 ( 



(1) 



where a is the semi-major axis of the orbit, i is the inclination, 
and 9 a is the eclipse semi-angle, we find that the observed A<j) = 
0.2 is inconsistent with a luminosity class III or V primary. We 
therefore conclude that indeed the primary is an 08.51 supergiant 
sta r (with mass, radius , luminosity, effective temperature given 
by iMartins etafl d2005l) . M OB = 31.5M , Rob = 21.4R , 
log Lob /Lq = 5.65, T efi = 32 274 K). Therefore, the distance 
to IGR J16418— 4532 is ~ 13 kpc dRahoui et alj|2008[) . 

During the eclipse of IGR J16418 — 4532, we detect a non-zero 
flux with the XRT (L x ~ 3 x 10 34 ergs _1 , see Fig. |4j. Detec- 
tions of HMXBs during the eclipse are quite common (see e.g. 
Ivan der Meer et alj|2005t IClark et alj[l98a ISako et al.|[l999l) and 
can be ascribed to X-ray photons scattered by the stellar wind of 
the OB star and discrete fluorescent and recombination lines from 
highly ionized gas around the compact object. The decrease in X- 
ray flux, when the compact object is approaching the eclipse, is 
related to the increase of the column density along the line of sight, 
through the stellar wind produced by the OB star. From obser- 
vations of other wind-fed X-ray pulsar s, such as 4U 1700—377 , 
Vela X-l, OAO 1657-415 (see e .g. Ivan der Meer et alj 120051 ; 
ISako et al.ll999] ; lAudlev et alj2006h the absorption column density 
during the eclipse ranges from ~ 10 23 cm -2 up to « 10 24 cm -2 . 

We can compare the X-ray luminosity and absorption column 
density of IGR J16418— 4532 during the eclipse with those cal- 
culated by Ivan der Meer et all J2005I) for 4U 1700-377, a wind- 
fed HMXBs with properties similar to IGR J16418-4532 (orbital 
period P or b = 3. 4 d, co mpanion star of spectral type 06.5Iaf). 
Ivan der Meer et al.l d2005h observed 4U 1700-377 with XMM- 
Newton during an eclipse; they found an X-ray luminosity of 
~ 7 x 10 34 ergs _1 (in agreement with the X-ray luminosity of 
IGR J16418 — 4532 during the eclipse) and a column density of 
^ 10 24 cm- 2 . 

The XRT light curve (Figs. [2] and [3} shows a large dynamic 
range, up to 1-2 orders of magnitude within a single snapshot and 
about 3 orders of magnitude overall (even when the low-intensity 
state is excluded). While this dynamic range is lower than that ob- 
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Figure 6. Observed cumulative luminosity distribution (solid black line) 
and relative uncertainties (dashed black lines), compared with a possible 
solution (red dot-dashed red line; see text). 



Table 3. Wind parameters values for IGR J 16418-4532. 



Parameter 


min. value 


max. value 


Mtot 


4 X 10" 7 M Q yr- 1 


8 X 10- 7 M Q yr- 1 


/ 


0.85 


0.95 


Voo 


800 kms -1 


1300 km s" 1 





0.8 


1.2 


c 


0.8 


1.6 


7 


-1.5 


-0.5 


Aia 


5 x 10 16 g 


5 x 10 19 g 


M b 


5 x 10 19 g 


10 21 g 



served in the classical dChatv et alj201Ch SFXTs, which reach 4-5 
orders of magnitude in dynamic range (e.g. IGR J17544— 2619), it 
is indeed typical of intermediate SFXTs (e.g., IGR J18483— 03 1 1). 
The 2011 February flare (Fig. [2] and dark grey filled circles at 
phase 0.27 in Fig. [4j, in particular, reached ~ 2 x 10 -9 erg 
cm~ 2 s _1 (2-10keV, unabsorbed), corresponding to a luminos- 
ity of ~ 4 x 10 37 ergs -1 , which is typical of SFXT outbursts as 
we have observed them wi th Swift/XKT (e.g., iRomano et alj200"8l 
l2009al ; ISidoli et al.ll2009bl l3). 

On the other hand, while the typical SFXT light curve is 
characterized by a faint flaring activity with luminosity of 1 33 — 
10 34 erg s -1 dSidoli et al.ll2008l ; lRomano et alj2009bll2011al) . spo- 
radically interrupted by flares with duration of a few hours and 
peak luminosity of 10 36 - 10 37 ergs -1 , IGR J16418— 4532 ex- 
hibits a flaring activity ranging from ~ 10 35 ergs -1 up to ~ 
4 x 10 37 erg s _1 (also see Walter et al. 2006 and the light grey curve 
in Fig. [4j . 

To reproduce the observed X-ray va riability, we applied the 
clumpy wind model of iDucci et alj ( 120091) . which we briefly sum- 
marize here. This model assumes that the wind of OB supergiants is 
inhomogeneous, composed of dense clumps surrounded by a hotter 
and homogeneous wind. Clumps follow a power-law mass distribu- 
tion 



-wo -*(£)' 



(2) 



where k = (1 - Q M~ ( /(Ml c - M a 1_c ) is the normaliza- 
tion constant, M c i is the mass of a clump, and M a and Mb de- 
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fine the extrema of the mass range. Assuming a spherical geometry 
for the clumps, the initial clump dimension distribution is given by 
Nm c1 oc H\ clumps s -1 , where R c \ is the radius of the clump. 
The clump mass loss rate M c \ is related to the total mass loss rate 
Mtot by means of / = M c \/Mtot, the fraction of wind mass in the 
form of clumps. Clumps and the homogeneous wind follow the /3- 
velocity law v(r) = Voo (l - 0.9983^a) /3 ^Castor et al.ll 19751) . 
where /3 is a constant in the range ~ 0.5 — 1.5, Rob is the radius 
of the supergiant and Voo is the terminal wind velocity. 

IGR J16418-4532 was intensely monitored by Swift/XKT for 
two orbital cycles (MJD 55755-55766). Each snapshot has a du- 
ration of ~ 1 ks, for a total of up to ~ 5 ks per day, which is 
the typical timesca le for the flare durations observed in SFXTs 
dSguera et al]|2006h . For this reason, we cannot establish the num- 
ber of observed flares and their duration, therefore we cannot use 
this information to compare the observ ed lightcurve with t hat cal- 
culated with the clumpy wind model of Du cci et al. j2009D. as has 
been done in prev ious work iDucci et alj2009l ; lRomano et alj20icl ; 
Ducc i et al .2010). However, it is possible to compare the observed 
cumulative luminosity distribution out of eclipse with that calcu- 
lated with the clumpy wind model. For the former, the Swift/XKT 
count rate, measured in the 0.3-10 keV, was converted to the 0.1- 
10 keV luminosity u sing the average spectra l parameters obtained 
bv lWalter etall j2006l) and lDucci et alj hold) , and assuming a dis- 
tance of 13kpc. For the latter, we adopted the physical parame- 
ters for the 08.51 primary and the neutron star secondary reported 
above, and t he force multip l ier pa rameters (see iDucci et alj|2009h 
obtained by IShimada eTail J 19941) for a 08.51 star (k = 0.375, 
q = 0.522, 8 = 0.099). 

Fig. [6] shows the Swift/XKT cumulative luminosity distribu- 
tion (black solid line) with its 90% uncertainties (black dashed 
lines). Fig. [6] also shows the cumulative luminosity distribution 
(red dot-dashed line) calculated with the clumpy wind model by 
adopting the following parameter values: mass loss rate Mtot = 
6 x 10~ 7 M o yr" 1 , terminal velocity = 1000 kms -1 , /3 = 1, 
fraction of mass lost in clumps / = 0.9, mass distribution power- 
law index £ = 1.1, power-law index of the initial clump dimension 
distribution 7 = — 1, minimum clump mass Af a = 5 x 10 18 g and 
maximum clump mass Mb = 10 20 g. Further acceptable solutions 
can be found assuming the wind parameters in the ranges reported 
in TablefJ] 

The wind parameters obtained for IGR J16418 — 4532 are 
roughly in agreement with those of HMXBs previously s tudied (see 
IDucci et al.ll2009l : lRomano et alJl20ld : lDucci et alfeOldl) . with the 
exception of / and Voo . The obtained low value of the terminal ve- 
locity can be explained as follows: when the wind is highly ionized 
by the X-ray photons emitted by the neutron star, the high ioniza- 
tion modifies the dynamics of the line-driven stellar wind of OB 
supergiants (especially in close binary systems), leading to a re- 
duc tion of the wind velocity in the direction of the neutron star (see 
e.g. IStevens & KallmarJ [l990l and references therein). The mass 
loss rate is lower than the typical mass loss ra te of Q8.5I stars , 
which is of the order of « 2 x 10 -6 M yr" 1 Jvink et al.ll2000h . 
but in agreement with the hypothesis that the mass loss rates de- 
rived from homogeneous-wind model measurements with the Ha 
method are overestimated by a f actor of 2-10 if the wind is clumpy 
(see e.g. iLepine & Moffatll2008h . However, we cannot exclude the 
possibility that the X-ray variability of IGR J16418-4532 is not to- 
tally due to the accretion of an inhomogeneous wind. In fact, other 
mechanisms could be at work, like e.g. transient accretion discs 
(Taametal. 198 8) and intermittent accr etion flow onto a neutron 
star fcamb et alJl977l : rDucci et aljioiol) . 



Our results, based on our intense monitoring of 
IGR J1641 8 — 4532, show that its X-ray properties (dynami- 
cal range, flaring activity and peak luminosity) are consistent with 
an intermediate SFXT nature. 
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